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Abstract. We present a weak lensing analysis of the double cluster system Abell 222 and Abell 223. The lensing reconstruction 
shows evidence for a possible dark matter filament connecting both clusters. The case for a filamentary connection between 
A 222/223 is supported by an analysis of the galaxy density and X-ray emission between the clusters. Using the results of 
A'-body simulations, we try to develop a criterion that separates this system into cluster and filament regions. The aim is to find 
a technique that allows the quantification of the significance of (weak lensing) filament candidates in close pairs of clusters. 
While this mostly fails, the aperture quadrupole statistics (Schneider & Bartelmann 1997) shows some promise in this area. 
The cluster masses determined from weak lensing in this system are considerably lower than those previously determined from 
spectroscopic and X-ray observations (Dietrich et al. 2002; Proust et al. 2000; David et al. 1999). Additionally, we report the 
serendipitous weak lensing detection of a previously unknown cluster in the field of this double cluster system. 
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1. The cosmic web 

The theory of cosmic structure formation predicts through N- 
body simulations that matter in the universe should be con- 
centrated along sheets and filaments and that clusters of galax- 
ies form where these intersect (e.g. Klypin & Shandarin 1983; 
Davis et al. 1985; Bertschinger & Gelb 1991; Bond et al. 
1996; Kauffmann et al. 1999). This filamentary structure, of- 
ten also dubbed "cosmic web", has been seen in galaxy red- 
shift surveys (e.g. Joeveer et al. 1978; de Lapparent et al. 1986; 
GiovanelH et al. 1986; Geller & Huchra 1989; Vogeley et al. 
1994; Shectman et al. 1996) and more recently by Baugh et al. 
(2004) and Doroshkevich et al. (2004) in the 2dF and SDSS 
surveys, and at higher redshift by MoUer & Fynbo (2001). 
Observational evidence for the cosmic web is recently also 
coming from X-ray observations. E.g. an X-ray filament be- 
tween two galaxy cluster was observed by Tittley & Henriksen 
(2001). Nicastro (2003) and Zappacosta et al. (2002) reported 
possible detections of warm-hot intergalactic medium fila- 
ments. 

Because of the greatly varying mass-to-light (M/L) ratios 
between rich clusters and groups of galaxies (Tully & Shaya 
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1999) it is problematic to convert the measured galaxy den- 
sities to mass densities without making further assumptions. 
Dynamical and X-ray measurements of the filament mass will 
not yield accurate values, as filamentary structures are probably 
not virialized. Weak gravitational lensing, which is based on 
the measurement of shape and orientation parameters of faint 
background galaxies (FBG), is a model-independent method 
to determine the surface mass density of clusters and fila- 
ments. Due to the finite ellipticities of the unlensed FBG ev- 
ery weak lensing mass reconstruction is unfortunately an in- 
herently noisy process, and the expected surface mass density 
of a typical filament is too low to be detected with current tele- 
scopes (Jain et al. 2000). 

Cosmic web theory also predicts that the surface mass den- 
sity of a filament increases towards a cluster (Bond et al. 1996). 
Filaments connecting neighboring clusters should have surface 
mass densities high enough to be detectable with weak lensing 
(Pogosyan et al. 1998). Such filaments may have been detected 
in several recent weak lensing studies. 

Kaiser et al. (1998) found a possible filament between two 
of the three cluster in the z - 0.42 supercluster MS 0302-1-17, 
but the detection remains somewhat uncertain because of a pos- 
sible foreground structure overlapping the filament and possi- 
ble edge effects due to the gap between two of the camera chips 
lying on the filament. Also, Gavazzi et al. (2004) could recently 
not confirm the presence of a filament in this system. Gray et al. 
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(2002, G02) claim to have found a filament extending between 
two of the three clusters of the Abell 901/902 supercluster, but 
the significance of this detection is low and subject to possible 
edge effects, as again the filament is on the gap between two 
chips of the camera. Clowe et al. (1998) reported the detection 
of a filament extending from a high-redshift (z = 0.809) cluster. 
Due to the small size of the image it is unknown whether this 
filament extends to a nearby cluster 

1.1. The A 222/223 system 

A 222/223 are two Abell clusters at z ^ 0.21 separated by 
~ 14' on the sky, or ~ 2800/z^g kpc, belonging to the Butcher 
et al. (1983) photometric sample. Both clusters are rich having 
Abell richness class 3 (Abell 1958). The Bautz-Morgan types 
of A 222 and A 223 are II-III and III, respectively. While these 
are optically selected clusters, they have been observed by 
ROSAT (Wang & Ulmer 1997; David et al. 1999) and are con- 
firmed to be massive clusters. A 223 shows clear sub-structure 
with two distinct peaks separated by ~ 4' in the galaxy distri- 
bution and X-ray emission. We will refer to these sub-clumps 
as A 223-S and A 223-N for the Southern and Northern clump, 
respectively. A 222 is a very elliptical cluster dominated by two 
bright elliptical galaxies of about the same magnitude. 

Proust et al. (2000) published a list of 53 spectra in the field 
of A 222/223, 4 of them in region between the clusters (here- 
after "intercluster region") and at the redshift of the clusters. 
Later Dietrich et al. (2002, D02) reported spectroscopy of 183 
objects in the cluster field, 153 being members of the clusters 
or at the cluster redshift in the intercluster region. Taking the 
data of Proust et al. (2000) and D02 together, 6 galaxies at the 
cluster redshift are known in the intercluster region, establish- 
ing this cluster system as a good candidate for a filamentary 
connection. 

1.2. Outline 

This paper is organized as follows. We describe observations of 
the A 222/223 system in Sect. 2. Our weak lensing analysis of 
this double cluster system is presented in Sect. 3; we compare 
this to the light (optical and X-ray) distribution in Sect. 4. We 
find possible evidence for a filamentary connection between 
the two clusters and try to develop a statistical measure for the 
significance of a weak lensing detection of a filament in Sect. 5. 
Our results are discussed in Sect. 6. 

Throughout this paper we assume a Q.^ = 0.7, On, = 
0.3, H() - IQIijo kms"' Mpc"' cosmology, unless otherwise 
indicated. We use standard lensing notation (Bartelmann & 
Schneider 2001) and assume that the mean redshift of the FBG 

isZpBG = 1- 

2. Observations of the A 222/223 system 

Imaging of the A 222/223 system was performed with the Wide 
Field Imager (WFI) at the ESO/MPG 2.2 m telescope. In total, 
twenty 600 s exposures were obtained in /?-band in October 
2001 centered on A 223, eleven 900 s /?-band exposures were 
taken in December 1999 centered on A 222. The images were 



taken with a dithering pattern filling the gaps between the chips 
in the co-added images of each field. 

The /?-band data used for the weak lensing analysis is sup- 
plemented with three 900 s exposures in the B- and V-band cen- 
tered on each cluster taken from November 1999 to December 
2000. The final B- and V-band images have some remaining 
gaps and regions that are covered by only one exposure and - 
due to the dithering pattern - do not cover exactly the same 
region as the /?-band images. 

The reduction of the /?-band image centered on A 222 is de- 
scribed in detail in D02. The /?-band image centered on A 223 
was reduced in the same way. The B- and V-band data was re- 
duced using the GaBoDS pipeline (Schirmer et al. 2003; Erben 
et al. 2005), using Astrometrix' with the USNO-A2 catalog 
(Monet et al. 1998) for the astrometric calibration and SWarp^ 
for the co-addition of the individual dithered images and chips. 
The B- and V-band pointings were co-added into a single frame 
for each color. The PSF properties of the /?-band pointings were 
so different that they were used separately for the lensing anal- 
ysis. The seeing of the co-added /?-band images is 0'.'9 and 0'.'8 
for the A 222 and A 223 pointings, respectively. 

The /?-band image centered on A 222 was photometri- 
cally calibrated using Landolt standard fields and corrected for 
galactic extinction (Schlegel et al. 1998), while the zero-point 
of the /?-band image centered on A 223 was fixed to match 
the magnitudes of objects in both fields. Because the B- and 
V-band data were known to be taken under non-photometric 
conditions, the red cluster sequence was identified in a color- 
magnitude diagram and its color adjusted to match those ex- 
pected of elliptical galaxies at the cluster redshift, using passive 
evolution and /T-correction on the synthetic galaxy spectra of 
Bruzual & Chariot (1993), to account for the additional atmo- 
spheric extinction. 

Due to the greatly varying coverage of the fields, it is diffi- 
cult to give a limiting magnitude for the co-added images. The 
number counts stop following a power law at 22.5-23.0 mag 
for the B- and V-band images and at 24 mag for the i?-band 
images. 

3. Lensing analysis 

3.1. Lensing catalog generation 

Starting from the initial SExtractor (Bertin & Arnouts 1996) 
catalog which contains all objects with at least 3 contiguous 
pixels 2a- above the background, we measured all quantities 
necessary to obtain shear estimates from the KSB (Kaiser et al. 
1995) algorithm. For this, we closely followed the procedure 
described in Erben et al. (2001). 

From the KSB catalog a catalog of background galaxies 
used for the weak lensing analysis was selected with the follow- 
ing criteria. Objects with signal-to-noise (SNR) < 2, Gaussian 
radius rg < 0'.'33 or rg > 1''19, or corrected ellipticity s > 0.8 
were deleted from the sample. Objects brighter than R < 22 
were rejected as probable foreground objects, while all objects 
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Fig. 1. Reliability tests of the shear estimates of objects observed in 
the overlapping region of the two R-band pointings. All objects with 
shear estimates are plotted as light dots, objects surviving our various 
selection criteria, detailed in the text, are plotted as crosses. Top left: 
Scatter plot of the ei component estimates from the exposures cen- 
tered on A 222 versus the one centered on A 223. The diagonal line is 
not a fit but only represents the ideal relation. Although we delete ob- 
jects with corrected ellipticity |£| > 0.8 from our final lensing catalog, 
some objects with |£| > 0.8 are marked with crosses in this plot. This is 
because the mean ellipticity of their two measurements, which we em- 
ploy in our selection and lensing analysis is below the chosen cut-off 
level. These objects are strongly down-weighted and their exclusion 
would not lead to significant differences in the lensing analysis. Top 
right: Same for the £2 component. Bottom left: Dependence of the ab- 
solute value of differences of the shear estimators |Ae| on the apparent 
magnitude of the object. As expected, fainter objects have less reli- 
able shear estimates. Bottom right: This panels shows the correlation 
between |Ae| and the weighting scheme we employed. Objects with 
higher weights have more reliable shear estimates. 



with R > 23 were kept as likely background galaxies. Objects 
between 22 < R < 23 with colors matching those of galaxies 
at redshift z < 0.5, -0.23 < (V - R) - 0.56 x{B-V) < 0.67, 
0.5 < B-V < 1.6 were not used for the lensing catalog. Objects 
not detected in the B-band image were kept if V - R > 1.0. 
The final catalog has 25940 galaxies, or 13.5 galaxies arcmin"^, 
without accounting for the area lost to masked reflection rings, 
diffraction spikes, and tidal tails. 

The large overlap between the two /?-band images allows 
us to test the reliability of the shear estimates and the validity 
of the weighting scheme we will employ in the lensing analy- 
sis. We perform these tests separately for the set of all objects 
found in the unmasked part of the overlap region, and the set of 
objects left after performing the various cuts described in the 
previous paragraph. The top panels of Fig. 1 show a compar- 
ison of the shear estimates of objects observed in the overlap 
region of the two /?-band pointings. Overall, the two indepen- 
dent shear estimates agree but show a broad scatter around the 



ideal relation. For the set of all galaxies, we find that the mean 
of the differences between the two measurements is -0.01 for 
the fii component and 0.00 for the 82 component. The standard 
deviation is 0.20 in each component. Erben et al. (2003) found 
an rms scatter of 0. 16 between two different lensing analyses of 
their data. Our value seems to indicate that the additional scat- 
ter introduced by the independent observations is small com- 
pared to the uncertainties intrinsic to the shear estimation pro- 
cedure. The bottom left panel of Fig. 1 shows the dependence 
of the absolute values of differences of the shear estimators |Ae| 
on the apparent/? -band magnitude. As one expects, the reliabil- 
ity of the shear estimates drops dramatically for fainter objects. 
Because these are the objects we keep in our lensing catalog the 
rms scatter between the two shear estimates increases to 0.25 
per component for the galaxies kept in our lensing catalog. The 
mean for the set of galaxies in our lensing catalog stays almost 
unchanged at 0.01 in both components, showing that, while the 
shear estimates become noisier, no systematic differences be- 
tween both images are present. 

In the lensing reconstruction and the aperture mass maps 
(Schneider 1996) we will assign a weight to each shear estima- 
tor. The weight is computed by 



w 



= «--D" 



(1) 



where o-g^^ is the intrinsic 2-d ellipticity dispersion and cr^ is 
the error estimate of the initial ellipticity measurement of the 
galaxy. We set (Xg^p = 0.38 which is typically found in ground- 
based weak lensing observations (T. Erben, private commu- 
nication; e.g. Clowe & Schneider (2001) who find a value of 
^EiT, - 0-42). cTg is computed from the uncertainty of the mea- 
surement of the quadrupole moment of the galaxy in the image. 
While both quantities are not independent - cr^^^ is of course 
increased by higher errors in the initial ellipticity measurement 
- their relation is very complex and not readily quantifiable 
in the KSB algorithm. As a consequence, galaxies with low 
CTg probably receive less weight than they should in an ideal 
weighting scheme. The large overlap and the high number of 
objects detected in both frames would enable us to study the 
shear estimation procedure in more detail and probably find a 
better weighting scheme than the one used in this work. This 
is, however, beyond the scope of this paper. 

The bottom right panel of Fig. 1 shows the correlation be- 
tween our weights (normalized to be < 1) and |Ae|. This verifies 
that galaxies with more reliable shear estimates have a higher 
weight in the generation of the lensing maps, although the large 
variations in |Ae| only correspond to small relative changes in 
the weight w. The shear estimates with the highest weight, not 
part of our lensing catalog, are those which we reject as proba- 
ble foreground objects because they are too bright. 

3.2. Weak lensing reconstruction 

Based on the lensing catalog described in the previous section, 
the weak lensing reconstruction in Fig. 2 was performed using 
the Seitz & Schneider (2001) algorithm adapted to the field 
geometry with a o" = 1!75 smoothing scale on a 214 x 200 
points grid. 
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Fig. 2. Weak lensing surface mass density contours overlaid on the i?-band mosaic observed with the Wide-Field imager at the ESO/MPG-2.2m 
telescope. The shear field was smoothed with a cr = 1^75 Gaussian, corresponding to the diameter of the circle at the lower left corner. Each 
contour represents an increase in k of 0.005 (~ 1.6 x 10'' Ayo Mq Mpc"^, assuming zfbg = 1) above the mean k at the edge of the field. 



Both clusters are well detected in the reconstruction, the 
two components of A 223 are clearly visible, and the ellipti- 
cal appearance of A 222 is present in the surface-mass density 
map. The strong mass peak West of A 223 is most likely asso- 
ciated with the reflection ring around the bright V - 7.98 mag 
star at that position. Although the prominent reflection ring was 
masked, diff'use stray light and other reflection features are vis- 
ible, extending beyond the masked region, well into A 223, 
probably being the cause of this mass peak. 

The peak positions in the weak lensing reconstruction are 
off-set from the brightest galaxy in A 222 and the two sub- 
clumps of A 223. The centroid of the mass of A 222 is 57" 
South-East of the brightest cluster galaxy (BCG); the mass cen- 



troids of A 223-N and A 223-S are 86" and 37" away from the 
BCGs of the respective sub-clumps. 

To estimate the significance of these off-sets we performed 
lensing simulations with singular isothermal sphere (SIS) mod- 
els of various velocity dispersions. The SIS models were put 
at the cluster redshift of z = 0.21; catalogs with a random 
distribution of background galaxies with a number density of 
15 arcmin"^, and 1-d ellipticity dispersion of cr^jp - 0.27 were 
created for 200 realizations. The shear of the SIS models was 
applied to the galaxy ellipticity of the catalog. Weak lensing re- 
constructions based on these catalogs were performed with the 
smoothing scale set cr = 1 !75 to match the smoothing of our 
real data. Due to the lower SNR for the 550 km s ' SIS, the sim- 



J. p. Dietrich et al.: Weak lensing study of dark matter filaments and application to the binary cluster A 222 and A 223 



0.8 



0.6 



-o 

V 




log-likelihood function (Schneider et al. 2000) 



y In p,(ei\gm) 



(2) 



0.4 



0.2 - 



40 60 80 100 

ceriLroid off-set/arcsec 

Fig. 3. Cumulative fraction of off-set of the reconstructed centroids 
from the real centroid. The curves display the probability of find- 
ing a reconstructed centroid of an SIS with a velocity dispersion 
of 550 km s"' (continuous), 700 km s"' (dotted), 850 kms"'(short 
dashed), 1000 kms"'(long dashed), and 1150 kms"'(dashed-dotted) 
from the real centroid position. The SIS was put at a redshift of 
z = 0.21; the number density of the input catalog was 15 arcmin"^. 



ulations yielded only 198 reliable centroid positions, while the 
centroid positions of the more massive SIS could be reliably 
determined in all 200 realizations. Fig. 3 shows the cumulative 
fraction of reconstructed peaks found within a given distance 
from the true centroid position. 

These simulations show that the observed off-set of A 223- 
S is compatible with the statistical noise properties of the re- 
construction. The off-set of A 222, using the the velocity de- 
termination of the SIS models fitted below, is significant at 
the 2-3cr level. The off-set of A 223-N cannot be explained 
with the statistical noise of the reconstruction alone. It is, how- 
ever, likely that the observed significant off-sets are not real 
but linked to the influence of the bright star and its reflection 
rings West of A 223. Although objects coinciding with this re- 
flection ring were excluded from the catalog, the presence of 
a strong mass peak on the position of the bright reflection ring 
is a clear indication that the shear estimates are affected by the 
weaker reflection features which are too numerous and large to 
be masked. It is difficult to guess how these reflections could 
contribute to the observed peak shifts. We found that varying 
the size of the masked region did affect the strength of the peak 
on the reflection ring but left the off-sets of the cluster peaks es- 
sentially unchanged. Still, it is noteworthy that the mass peaks 
are shifted preferentially away from the star. 

To avoid the mass-sheet degeneracy we estimate the clus- 
ter masses from fits of parameterized models to the shear cat- 
alog. The fits were performed minimizing the negative shear 



over the Ny galaxy images to obtain the parameter set most con- 
sistent with the probability distribution PsisilgiOd) of lensed 
galaxy ellipticities. See Schneider et al. (2000) for a detailed 
discussion of this maximum likelihood method for parameter- 
ized models. We fitted more than one mass profile simultane- 
ously. Compared to a single model fit, this reduces the influ- 
ence of the other cluster on the fitting procedure and result. 
Galaxies within distances 6 < 6^i„ = 3' from the centers of 
the models were ignored when fitting SIS models. Assuming a 
typical SIS this corresponds to roughly 10 Einstein radii and is 
enough to assume that all galaxies used in the fitting procedure 
are in the weak lensing regime. Ignoring galaxies close to the 
cluster centers also reduces the contamination with faint clus- 
ter galaxies. As a first approach we fit two SIS, one centered 
on the BCG of A 222, the other centered on the line connecting 
the BCGs of the two sub-clumps of A 223. The best-fit mod- 
els in this case have velocity dispersions of 716^^4 km s"' and 
804^^^ km s"' , respectively. This is considerably lower than the 
spectroscopic velocity dispersions of D02 of 1014^^'| km s^'for 
A 222 and 1032^^^ km s^'for A 223. It is also lower than the ve- 
locity dispersions derived from X-ray luminosities (David et al. 
1999) and the Lx - cr relation of Wu et al. (1999), which are 
845 - 887 kms-'and 828 - 871 km s"', respectively (D02), but 
the value for A 223 is consistent within the error bars. The er- 
ror bars of the individual velocity dispersion were computed 
from 2A/, where the velocity dispersion of one component was 
kept fixed at its best-fit value to give the eiTors estimate for the 
other component. The two component fit has a significance of 
7.9cr over a model without mass. Joint confidence contours are 
displayed in Fig. 4. A three component model with an SIS cen- 
tered on each BCG has a lower significance over a zero mass 
model than the two SIS model and does not fit the data better. 

Because both clusters are elliptical and the masses deter- 
mined from SIS fits differ strongly from those derived by D02, 
one might assume that fitting elliptical mass profiles yields a 
more accurate estimate of the cluster mass. To test this, we fit- 
ted singular isothermal ellipse models to the clusters. It turned 
out that the 6 parameter fit necessary to model both cluster si- 
multaneously was very poorly constrained and the fit proce- 
dure was not able to reproduce the orientation of the clusters. 
Results strongly depended on the initial values chosen for the 
minimization routines. 

The best-fit NFW models have raoo = 1276!}^^/i7,J kpc, c = 
3.4 and rzoo = 1546^\f^h-^ kpc, c = 1.2 for A 222 and A 223, 
respectively, excluding shear information at distances < 1 ^5 
from the cluster centers. The NFW models have a significance 
of 5.2cr over a model with no mass. Fig. 5 shows confidence 
contours for the NFW fits to the individual clusters, computed 
from 2A/, while keeping the best-fit parameters for the other 
cluster fixed. We summarize the derived cluster properties in 
Tab. 1. 

As we see from the left panel in Fig. 5, it can be difficult 
to obtain reliable concentration parameters from weak lensing 
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Fig. 5. Confidence contours for the best-fit NFW parameters for A 222 {left panel) and A 223 {right panel). The contours are drawn at the same 
levels as in Fig. 4. 
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Fig. 4. Combined confidence contours for the SIS velocity disper- 
sion of A 222 and A 223. The contour lines are drawn for 2AZ = 
(2.3,4.61,6.17,9.21,11.8,18.4), corresponding to the 63.8%, 90%, 
95.4%, 99%, 99.73%, and 99.99% confidence levels, respectively, un- 
der the assumption that the statistics is approximately Gaussian. 



data. The reason is that the shear signal is mostly governed by 
the total mass inside a radius around the mass center. Only in 



Table 1. Summary of the cluster properties derived from spectroscopy. 
X-ray, and weak lensing observations. 
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the cluster center the shear profile carries significant informa- 
tion about the concentration parameter. For example if we set 
^min = 3' - like we did for fitting SIS models - in the minimiza- 
tion procedure, M200 remains essentially unchanged while the 
concentration factor can increase dramatically. The best fit pa- 
rameters for A 222 in this case are r2oo - 1238/;^^ kpc, c - 7.8. 
If we choose the radius fl^in inside which we ignore galaxies 
too big, typical values for the scale radius rs - r2oo/c are con- 
tained within this radius, c is then essentially unconstrained, 
i.e. the minimization procedure cannot anymore distinguish be- 
tween a normal cluster profile and a point mass of essentially 
the same mass. 

The situation is different for A 223. The two sub-clumps 
are separated by ~ 4'. This means that even ignoring shear in- 
formation within the larger 0min = 3' radius, the outer slopes 
of the sub-clumps are outside 0n,in and the determination of the 
concentration parameter gives a tight upper bound and does not 
change as dramatically when the minimization is performed 
only with galaxies further away from the cluster center as is 
the case in A 222. Because the shear outside O^i^ is effectively 
that of an averaged mass profile inside ^min the measured con- 
centration parameter is very low. 

The projected cluster separation is marginally smaller than 
the sum of the virial radii (rvir ~ ^200) derived from the shear 
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analysis; r2oo(A 222) + r2oo(A 223) = 2822+[^^/z7q' kpc. We 
have to emphasize that this is only the projected separation. 
D02 found redshifts of z = 0.2126 + 0.0008 for A 222 and 
z = 0.2079 ± 0.0008 for A 223. Assuming that both clus- 
ters participate in the Hubble flow with no peculiar velocity, 
this redshift difference of Az = 0.005 + 0.001 translates to a 
physical separation along the line of sight of (15 ± 'i)h^^ Mpc. 
Presumably, part of the observed redshift difference is due to 
peculiar velocities. In any case, it is more likely that the two 
clusters are physically separated and the virial radii do not over- 
lap. 

3.3. A possible dark matter filament 

Also visible in the reconstruction is a bridge in the surface 
mass density extending between A 222 and A 223. Although 
the signal of this possible filamentary connection between the 
clusters is very low, the feature is quite robust when the selec- 
tion criteria of the catalog are varied and it never disappears. 
Variations on the selection criteria of the catalog let the fila- 
ment shift a few arcminutes in the East- West direction. The 
filament strength also changes but on closer inspection this can 
be attributed almost entirely to variations in the mass-sheet de- 
generacy, which is fixed by setting the mean k at the edge of 
the field to zero. Although the field is big enough to assume 
that the clusters have no considerable contribution to the sur- 
face mass density at the edge of the field, this is a region where 
the K-map is dominated by noise. Small changes in the selec- 
tion criteria can change the value of the mass-sheet degeneracy 
by as much as kq - 0.02. This illustrates that a surface mass 
density reconstruction is not suitable to assess the significance 
of structures as weak as filaments expected from A^-body simu- 
lations. We try to develop methods to quantify the significance 
of this signal in section 5. 



3.4. Other mass peaks 

In addition to the cluster peaks several other structures are seen 
in the reconstruction in Fig. 2. Using the aperture mass statis- 
tics (Schneider 1996) with a 6.'4 filter scale we find that the 
peak ~ 13' SE of A 222 has a SNR of 3.5. This peak corre- 
sponds to a visually identified overdensity of galaxies. Fig. 7 
shows a V - R vs. R color-magnitude diagram of all galaxies 
in a box with 170' side length around the brightest galaxy in 
this overdensity. A possible red-cluster sequence (RCS) can be 
seen centered around V - R = 1.1, which would put this mass 
concentration at a redshift of z ~ 0.4. However, the locus of 
the RCS is so poorly defined that this estimate has a consider- 
able uncertainty. Assuming this redshift, the best-fit SIS model 
has a velocity dispersion of 728^J2okm s^'and a significance of 
3.2cr over a model without mass. The best-fit NEW model has 
'"200 = 1322/z^Q kpc and c = 3.3 and a significance of 2.8cr over 
a model with r2oo - 0, determined from Sx^- 

Fig. 8 shows a comparison of surface mass and luminos- 
ity density for this peak. Galaxies with 1.0 < V - R < 1.2 
were selected to match the tentative RCS from Fig. 7. The 
figure shows excellent agreement between the mass and light 
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Fig. 6. Shown above are SNR contours of the aperture mass statistics 
with a 6f 1 filter radius overlaid on the WFI i?-band image. The lowest 
contour is at 2.0, higher contours rise in steps 0.5. The mass peak ~ 13' 
SE of A 222 has a peak SNR of 3.5. The circle segment in the lower 
left corner has the same radius as the filter function. 
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Fig. 7. Color magnitude diagram of objects around the mass peak SE 
of A 222. A possible red-cluster sequence can be seen around V-R ~ 
1.1. 



contours, unambiguously confirming that this is a weak lens- 
ing detection of a previously unknown cluster. The off-set be- 
tween the mass centroid and the BCG, which is located at 
CK=01:38:12.1, ^=-13:06:38.2, is 27" and not significant for 
an SIS with a velocity dispersion of ~ 730 kms 'at a redshift 
ofz~0.4. 
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Fig. 8. Mass and light contours in the peak SE of A 222. Black solid 
lines are contours of the mass reconstruction in Fig. 2, white dashed 
lines denote the luminosity density of galaxies with l.O < V-R < 1.2. 



The mass concentration in the Western part of the possible 
filament reaches a peak SNR of 3.6 at a filter scale of 3'2. We 
do not find an overdensity in the number or luminosity density 
of galaxies at this position. None of the other mass peaks seen 
in the reconstruction, with exception of the one on the reflection 
ring, is significant in filter scales > 2!4. 

4. Comparison of mass and light in A 222/223 

Fig. 9 shows the number density and luminosity distribution in 
the A 222/223 system. The left panel shows the number density 
distribution of the color-selected (0.78 < V - R < 0.98) early- 
type galaxies. The contours lines indicate the SNR determined 
from bootstrap resampling the selected galaxies. It is evident 
that a highly significant overdensity of early-type galaxies ex- 
ists in the intercluster region. The right panel shows a compar- 
ison between luminosity (background gray-scale image) and 
surface mass density. In general, there is good agreement be- 
tween the two. We again note the off-sets of the mass centroids 
from the light distribution, which we attribute to the systemat- 
ics induced by the reflection ring. The elongation of A 222 is 
nicely reproduced in the reconstruction. A 222 is the dominant 
cluster in the luminosity density map, while in the mass recon- 
struction A 223 appears to be more massive. We note, however, 
that many of the bright E/SO galaxies in A 223 escape our color 
selection because they are bluer than expected for early-type 
galaxies at this redshift. This indicates a high amount of merger 
activity in this irregular system and most likely still collapsing 
system. As we fixed the colors such that the RCS matches the 
expected colors of early-type galaxies at the cluster redshift, 
this clearly shows that the bright central galaxies in A 223 are 
bluer than expected. 



The overdensity in number and luminosity density is not 
aligned with the dark matter filament candidate. We should, 
however, not forget that the position of the filamentary structure 
is somewhat variable with varying cuts to the lensing catalog. 
If the reflection features West of A 223 are indeed responsible 
for shifting the centroid positions of the massive clusters, their 
effect may be even stronger on such weak features as the mass 
bridge seen in the reconstruction. 

We estimate the cluster luminosities by measuring the R- 
band luminosity density of all galaxies within r2oo - as deter- 
mined from the lens models in Sect. 3.2 - in excess of the lumi- 
nosity density in a circle with 5' radius centered on (01:36:45.8, 
-13:07:25), which is an empty region in the SW of our field. 
A 222 has a luminosity Lr, ^,00 = (2-7 ± 0.4) x 10^^ h'^ Lq and 
A 223 has a luminosity of Lg^^.^^ = (5.6 + 0.8) x lO'^ h'^ Lq. 
Using the mass determined from the NFW profiles, this im- 
pHes rather low M/L ratios; M/Lr - lll^Uhjo Mq/Lq for 
A 222 and M/Lr = 95^11 ho Mq/Lq for A 223. The mass- 
to-light ratios increase to M/Lr - ll^'^'^g IijoMq/Lq and 
M/Lr = 131^11 ho Mq/Lq for A 222 and A 223, respectively, 
if the mass estimates from the SIS models within ^200 are used. 

The X-ray satellite ROSAT observed the pair of galaxy clus- 
ters on 16. January 1992 using the position sensitive propor- 
tional counter (PSPC). We extracted these data from the public 
ROSAT archive in Munich and analyzed the total integration 
time of 6780 seconds using the EXSAS software (Zimmermann 
et al. 1998). To avoid any confusion with diffuse soft X-ray 
emission and associated photoelectric absorption towards the 
area of interest, we focused our scientific interest on the upper 
energy limit of the PSPC detector. Using the pulse height in- 
variant channels 51 - 201 (corresponding to 0.5 keV < E < 
2. 1 keV) we calculated the photon image and the correspond- 
ing "exposure-map" according to the standard data reduction. 
We performed a "local" and a "map" source detection which in 
total yielded 42 X-ray sources above a significance threshold 
of ten. 

By visual inspection, we selected some X-ray sources lo- 
cated close to the diffuse X-ray emission of the intra-cluster 
gas and subtracted their contribution using the EXSAS task 
create/bg_image. This task subtracts the X-ray photons of 
the point sources and approximates the background intensities 
via a bi-cubic spline interpolation. Finally the X-ray data were 
smoothed to an angular resolution of 1 f75 using a Gaussian 
smoothing kernel. 

Contours for this final image are shown in Fig. 10. Detected 
X-ray sources kept in the final image are marked with circles; 
the subtracted unresolved sources are denoted by stars. The 
lowest contour line is at the 3cr level. Higher contours increase 
in steps of 2cr. Both cluster are very well visible. As already 
noted by Wang & Ulmer (1997), A 223-S is by far the domi- 
nant sub-clump in A 223 in X-ray. 

A bridge in X-ray emission connecting both clusters is seen 
at the 5cr level in this image. This possible filament is aligned 
with the overdensity of the number density of color selected 
galaxies but not with the filament candidate seen in the weak 
lensing reconstruction. 

The Eastern spur in the X-ray emission of A 223 is caused 
by a point source whose removal would cut significantly into 
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Fig. 9. Smoothed distribution of the number density (left panel) and the luminosity density (right panel). The smoothing was done with a 
cr = if 75 Gaussian to match the smoothing of the weak lensing reconstruction. The diameter of the circles at the lower left corners corresponds 
to the FWHM of the Gaussian. The contour lines in the left panel are significance contours starting at 5cr and rising in steps of Icr, the contour 
lines in the right panel are the surface mass density contours of the reconstruction. 



the cluster signal. We therefore decided to keep this source. 
Removing it does not influence the signal in the intercluster re- 
gion. The Northern extension of A 223 in the k map is blinded 
by the support structure of the PSPC window in the X-ray ex- 
posure. 

5. Quantifying filaments 

In Sect. 3 we found a possible dark matter filament extending 
between the clusters A 222 and A 223 but its reality and sig- 
nificance are not immediately obvious. In this section we try to 
develop a method to quantify the significance of weak lensing 
filament candidate detections. 

To quantify the presence of a filament and the significance 
of its detection, two problems must be solved. First the funda- 
mental question "What is a filament?" must be answered. How, 
for instance, is it possible to discriminate between overlapping 
halos of two galaxy clusters and a filament between two clus- 
ters? While in the case of large separations this may be easy 
to answer intuitively, it becomes considerably more difficult if 
the cluster separation is comparable to the size of the clusters 
themselves (see e.g. the left panel of Fig. 9). 

The second problem - quantifying the significance of a fil- 
ament detection - is rooted in the weak lensing technique. To 
avoid infinite noise in the reconstruction, the shear field must 
be smoothed (Kaiser & Squires 1993). This leads to a strong 
spatial correlation of the nns error in the reconstructed k map, 
making it difficult to interpret the error bars at any given point. 
Randomizing the orientation of the faint background galaxies 
while keeping their ellipticity moduli constant and performing 
a reconstruction on the randomized catalog allows one to assess 
the overall noise level as (f^^^^). While this can be used to de- 



termine the noise level, the mass-sheet degeneracy (Schneider 
& Seitz 1995) allows us to arbitrarily scale the signal. 

Statistics like the aperture mass (Schneider 1996) and aper- 
ture multipole moments (Schneider & Bartelmann 1997) allow 
the calculation of signal-to-noise rations for limited spatial re- 
gions and are thus well suited to quantify the presence of a 
structure in that region. Hence, to quantify the presence of a 
structure between two galaxy clusters, the aperture has to be 
chosen such that it avoids the clusters and is limited to the fil- 
ament candidate. This is of course closely related to the first 
problem. We will discuss in the following sections how aper- 
ture statistics could be used to determine the SNR of a possible 
dark matter filament. 

We use A^-body simulations of close pairs of galaxy clusters 
to find solutions to these problems. 

5.1. N-body simulations 

Since we are interested in developing a method for a very par- 
ticular mass configuration, it is desirable to work with A^-body 
simulations that could mimic as closely as possible the A 222 
and A 223 cluster system. This goal can be achieved by per- 
forming constrained A^-body simulations. 

Constrained realizations were first explored by Bertschin- 
ger (1987), and later presented in an elegant and simple formal- 
ism by Hoffman & Ribak (1991). Here we follow the approach 
of van de Weygaert & Bertschinger (1996) of the so-called 
Hoffman-Ribak algorithm for constrained field realizations of 
Gaussian fields. With this approach, the constructed field obeys 
the imposed constraints and replaces the unconstrained field. 

The cluster-bridge-cluster system intended to simulate re- 
sembles a quadrupolar matter distribution. It is therefore ex- 
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Fig. lO.To the left is an overlay of X-ray 
contours (white lines) over the WFl images 
of our field. The contour lines start at 3cr 
and increase in steps of 2cr. All detected X- 
ray sources are marked; circles correspond 
to sources kept, while the point sources ex- 
cised from the image are marked with stars. 
The X-ray image was smoothed with a cr = 
1175 Gaussian, corresponding to the diam- 
eter of the circle at the lower left comer. 
The black contours in the background are 
the significance contours of the number den- 
sity from the left panel of Fig. 9. 



pected that primeval tidal shear plays an important role in shap- 
ing such matter configuration (see van de Weygaert 2002 for a 
review), which must have been induced by tiny matter density 
fluctuations in the primordial universe. 

In moulding the observational data we have considered 
each one of the observational aspects and cosmological charac- 
teristics of the system. We have put constraints onto the initial 
constrained field to create the two clusters and the bridge in the 
following way: 

Two initial cluster seeds were sowed at the center of the 
simulation box, separated by a distance d. We imposed con- 
straints over the clusters themselves like peak height (1 con- 
straint), shape (3 constraints), orientation (3 constraints), pecu- 
liar velocity (3 constraints) and tidal shear (5 constraints). 

We have performed a set of 10 realizations in a periodic 
5O/I7Q Mpc box, with different combinations of the mentioned 
18 constraints. In all simulations we set hjo - 0.7, i.e. //q = 
50 km s '. All constraints were imposed over a cubic grid 
of 64 grid-cells per dimension. In all 10 realizations, all con- 
straints were defined on a Gaussian scale r^ of 4 h^o Mpc for 
both clusters. Because we are dealing with rich clusters, we 
have imposed a peak height fc = 3cro, where cro is the vari- 
ance of the smoothed density field (cro(r^) = {.fgfgY^^)- The 
other constraints considered were: oblate clusters with axis ra- 
tios A2IA1 - 0.9 and At,IAi - 0.8 and both major axes aligned 
with each other We have imposed a "weak" primordial tidal 
field in order to produce a realistic field around the clusters. 
The stretching mode of the tidal field was aligned along the 
same direction given by the major cluster axes. The compres- 
sional mode was set perpendicular to the bridge axis. This com- 



bination of constraints proved to be the most successful one in 
reproducing (in the linear regime) the configuration presented 
by the two Abell clusters. 

The initial particle displacements and peculiar velocities 
were assigned according to the Zel'dovich (1970) approxima- 
tion from the constrained initial Gaussian density field. The 
evolution of the linear constrained density field into the non- 
linear regime was performed by means of a standard P^M code 
(Bertschinger & Gelb 1991). The number of grid-cells used to 
evaluate the particle-mesh force was 128-', with a particle mass 
resolution of 3.3 x 1O'°M0. We selected 15 time outputs in or- 
der to follow the simulation through the non-linear regime, with 
a time output at redshift z - 0.21, to match the observed clus- 
ter redshift. Fig. 1 1 shows the most successful cluster-bridge- 
cluster configuration. 

To estimate the underlying smooth mass distribution from 
the result of the simulations, the particle distribution was 
smoothed using the adaptive kernel density estimate described 
by Pisani (1996, 1993). A comparison of the simulated and 
smoothed mass distribution can be found in Figs. 11 and 12. 
The surface mass density of all simulations was linearly scaled 

such that A-jnax - 1 ■ 

We have also computed the density field by means of 
the Delaunay Tessellation Field Estimator (Schaap & van de 
Weygaert 2000), which in principle offers higher spatial res- 
olution at both, dense and underdense regions in comparison 
with other fixed-grid or adaptive kernels density reconstruction 
procedures. 

The reconstructed DTFE surface mass density maps of the 
cluster-bridge-cluster system agree with those from the adap- 
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Fig. 11. Zoom in on the central 10 x 10 Mpc-//i^g of an A'-body simu- 
lation. Displayed is the projection of a slice of 2.5 Mpc//!7o thickness. 
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Fig. 12. Smooth density distribution of the data in the left panel 
from the adaptive kernel density estimate. The contours are at /c = 
{0.03,0.05,0.1,0.51. 




Fig. 13. Reconstruction of the mass distribution in Fig. 12 on a 
206 X 206 points grid. The scale of the axes is given in arcminutes. 
The contours mark an increase of k in steps of 0.025 above the mean 
of the edge of the field. The diameter of the circle in the lower left is 
equal to the FWHM of the Gaussian smoothing of the shear field. 



package.^ kappa2 stuff solves the Poisson equation 

vV(6') = 2K{e) 



(3) 



in Fourier space on the grid by means of a Fast Fourier 
Transformation (FFT) and returns - among other quantities - 
the complex shear y{ff) and the magnification /i(0). 

For the lensing simulation, catalogs of background galaxies 
were produced. Galaxies were randomly placed within a pre- 
defined area until the specified number density was reached. 
To each galaxy an intrinsic ellipticity was assigned from two 
Gaussian random deviates. Unless noted otherwise all simula- 
tions have 30 galaxies/arcmin^ and a one-dimensional elliptic- 
ity dispersion of a-g^^ - 0.2. 

To test the validity of our lensing simulation we performed 
a mass reconstruction of the simulation in Fig. 12 using the 
algorithm of Seitz & Schneider (2001). The smoothing scale 
in this reconstruction was set to 1.'3. The result is shown in 
Fig. 13. We see that the reconstruction successfully recovers 
the main properties of the density distribution; both clusters 
are clearly detected, their ellipticity and orientation agrees with 
that of the smoothed density field. A "filament" resembling that 
in Fig. 12 is also seen. We now have to find a way to determine 
the significance of its detection. 



tive kernel since in high-density regions, both methods give 
similar density estimates (Pelupessy et al. 2003). 

The lensing properties of the smoothed mass distribu- 
tion were computed on a 2048 x 2048 points grid using the 
kappa2 stuff program from Nick Kaiser's IMCAT software 



5.2. Fitting elliptical profiles to galaxy clusters 

In a first attempt to quantify filaments, we try to fit the galaxy 
clusters by elliptical mass profiles. We then define the filament 



http://www.ifa.hawaii.edu/~kaiser/imcat/ 
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as the part of the mass distribution which is in excess of the 
mass fitted by the ellipses. 

The profile we used is the non-singular isothermal ellipse 
of Keeton & Kochanek (1998). This profile has four parameters 
per cluster that need to be fitted: 

- The axis ratio 

- The core radius 

- The Einstein radius of the corresponding singular isother- 
mal sphere 

- The orientation of the ellipse 

The central position was fixed and taken to be the peak po- 
sition of the mass reconstruction. As in Sect. 3.2 we use the 
shear-log-likelihood function to fit more than one mass profile 
simultaneously. 

Various methods for multidimensional minimization are 
available. All programs used for fitting either used the Down- 
hill Simplex or Powell's Direction Set algorithms discussed in 
detail in Press et al. (1992). We could not find any systematic 
differences between the results of the two methods. In general, 
their results agreed quite well if the same initial values were 
used. 

Fig. 14 shows a fit of two non-singular isothermal ellipse 
(NIE) profiles to the simulation in Fig. 12. We see that the 
cluster are roughly fitted by the NIE profiles. Like in the case 
of fitting SIE models to A 222/223, the ellipticity of the origi- 
nal cluster is only poorly reproduced. Varying the initial values 
of the minimization procedure, gives comparable best-fit val- 
ues for the Einstein and core radius. The axis ratio and orien- 
tation of the ellipse are so strongly affected by the choice of 
initial values, that they have a profound impact on the surface 
mass density in the intercluster region. Especially, the orienta- 
tion is only poorly constrained. Generally, the fits overestimate 
the surface mass density in the intercluster region, fitting the 
"filament" completely away. The behavior of the fits to the sim- 
ulated data confirms our experience with fitting SIE profiles in 
the A 222/223 system. Letting the slope of the density profile 
vary does not remedy the problem. The shear log-likelihood 
function is rather sensitive to the slope of the density profile, 
but the ellipticity of the clusters is still poorly constrained. 

5.3. Using aperture multipole moments to quantify the 
presence of a filament 

Aperture multipole moments (AMM) quantify the weighted 
surface mass density distribution in a circular aperture. If it is 
possible to find a characteristic mass distribution for filaments 
and express it in terms of multipole moments, AMM can be 
used to quantify the presence of a filament. 

Fig . 1 5 illustrates with the help of a simple toy model of two 
galaxy clusters connected by a filament why one expects to find 
a quadrupole moment in an aperture centered on the filament. 
Fig. 16 illustrates that it is crucial not to choose the aperture 
too large. If the aperture also covers the clusters, a quadrupole 
moment will be measured even if no filament is present. 




Fig. 15. Simple toy model of two galaxy clusters connected by a fila- 
ment. A quadrupole moment is present in the aperture centered on the 
filament. 




Fig. 16. Toy model of two galaxy clusters without a filament, illustrat- 
ing why it is important to choose the correct size of the aperture. 



We choose a weight function 

U(9) 



lO otherwise , 



(4) 



to compute the aperture quadrupole moment as defined in 
Schneider & Bartelmann (1997, see also Appendix B). While 
this weight function is clearly not ideal as it does not closely 
follow the mass profiles of the simulated data, it is sufficient to 
identify all relevant features in quadrupole moment I2*^'l maps. 
Fig. 17 shows such \Q.^^^\ maps for the simulation in Fig. 12. 
In the quadrupole maps 0niax increases from 2'to 5'. The maps 
were computed on 55 x 55 points grid, so that each grid point 
is 1' X 1' big. Overlayed are the contours of the surface mass 
density of the reconstruction of Fig. 13. 

One clearly sees that the quadrupole moment between the 
clusters increases as the size of the apertures increases. This 
is of course expected and due to the growing portion of the 
clusters in the aperture, so that their large surface mass den- 
sity dominates the mass distribution. Fig. 17 also illustrates the 
problem of separating clusters from a filament. The virial radii 
of the clusters extend far beyond the mass contours of the clus- 
ters in most directions and thus beyond what can be detected 
with weak lensing. Due to the ellipticity of the clusters it is not 
obvious whether the projected mass extending out to the virial 
radius is part of the cluster or belongs to a filament. While we 
certainly can define that projected mass outside the virial radii 
belongs to a filament, the case is not clear for mass inside r2oo- 
The surface mass density contours in Fig. 17 seem to suggest 
that a signature of a filament is present and observable with 
weak lensing inside the virial radius. Because weak lensing 
only has a chance to detect filaments in close pairs of clusters 
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Fig. 14. Left: Fit of two NIE profiles to the simulation in Fig. 12. The contour lines increase from a: = 0.025 lo k = 0.25 in steps of 0.025. Right: 
Difference image of Fig. 12 and the left panel. The contours are at the same level as in the left panel. The filament was completely subtracted. 
The surface mass density in the intercluster region is mostly negative. 



whose separation is comparable to the sum of their virial radii, 
understanding this signature is important. 

In this context the two maps in the top panel which show 
the smallest overlap of the aperture with the virial radii are the 
most interesting. Noteworthy in Fig. 17 is also that the top pan- 
els show a quadrupole moment on a ring-like structure around 
each cluster center. This is indeed to be expected for all galaxy 
clusters because there is a non-vanishing quadrupole moment 
if the aperture is not centered on the cluster center, but some- 
where on the slope of the mass distribution. This now raises the 
question how we should distinguish the quadrupole moment 
present around any cluster from that caused by a filamentary 
structure between the clusters. 

To better understand the features visible in the 12*^^'! maps 
of the A^-body simulations we qualitatively examine the struc- 
tures of a quadrupole map for a system of two isolated clus- 
ters and two clusters connected by a filament using simple toy 
models. Fig. 18 shows noise-free quadrupole maps of two trun- 
cated NFW halos (Takada & Jain 2003) without (top panel) and 
with (bottom panel) a connecting filament for the same aperture 
sizes as in the top panel of Fig. 17. The halo centers and virial 
radii were chosen to match those in the A^-body simulation. To 
describe the filament we choose a coordinate system such that 
the X-axis runs through the halo centers and has its origin at the 
center between the clusters separated by a distance d and define 
the following function: 



Ka(x,y) = ko — — ^ — ; ^ 

{ylhY + 1 



d 



(5) 



where A:i < and ^2 chosen such that the maxima of the fourth 
order polynomial coincide with the halo centers. H{x) is the 
Heaviside step-function. 



As in the \Qf^^\ maps of the A^-body simulations a qua- 
drupole moment related to the slope of the clusters is visible. 
Already the smallest aperture (left panel) overlaps the truncated 
NFW halos and leads to a strong quadrupole moment in the in- 
tercluster region. This quadrupole moment is, however, much 
weaker than it is in the presence of a filament (bottom panel 
of Fig. 18). Additionally, the presence of a filament is indi- 
cated by a ring structure on which the quadrupole moment is 
lower This structure becomes more prominent with increas- 
ing filter radius. All aperture statistics act as bandpass filters 
on structure comparable in size to the filter radius. As the ring 
has a radius of ~ 7' it is better visible in the map generated 
from the larger filter. This structure is present only in the halo- 
filament-halo system and not in the halo-halo system, even for 
filter scales larger than those depicted in Fig. 18. This structure 
is also visible in the quadrupole maps of the A^-body simula- 
tions. It is well visible in the top right panel of Fig. 17 and less 
well visible but still present in the top left panel of Fig. 17. 
Thus, the quadrupole maps clearly indicate that the measured 
quadrupoles on the filament are not caused by a situation with- 
out filament like that illustrated in Fig. 16. Unfortunately, this 
ring structure, like the filament itself, is a visual impression 
that does not lend itself easily to a quantitative assessment of 
the presence of a filament. 



These toy models and the A'-body simulation show that the 
aperture quadrupole statistics is indeed sensitive to filamen- 
tary structures. The quadrupole moment of a halo-filament-halo 
system exceeds that of a pure halo-halo system. As the mea- 
sured quadrupole moment strongly depends on the size of the 
aperture, choosing the appropriate aperture is important. A de- 
composition of a halo-filament-halo system into halo and fil- 
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Fig. 17. Quadrupole moment of the simulation in Fig. 12. Overlayed 
are the contours of the mass reconstruction in Fig. 13 and large circles 
with radius r2m as determined from the 3-dimensional simulated data. 
1 2*^' I was computed in an aperture with radius Smax = 2' (top left), 
Omax - 3' (top right), &max = 4' (bottom left), 0„ax - 5' (bottom right). 
The circles at the lower left comers have the same radii as the filter 
function in the respective panel. The cluster on the left has a mass of 
M200 = 2.0 X lO'^M©, the one on the right of M200 = 1-4 x lO^M©. 



ament components could be provide an objective criterion for 
chosing the filter scale. 

While in the (failed) attempt to separate the clusters and 
the filament by fitting elliptical profiles to the clusters, the fila- 
ment was naturally defined as the surface mass density excess 
above the clusters, there is no criterion in the AMM statistics 
that defines cluster and filament regions. We try to develop such 
a criterion in the following section. 

5.4. Defining cluster and filament regions 

Fig. 19 shows a simple one-dimensional toy model of the 
mass distribution of a cluster with a filamentary extension. 
The model consists of the following components: We assume 
a cluster with a King profile. This is the solid line in Fig. 19. 
In all simulations we see that the clusters are not spherical but 
triaxial with their major axes oriented approximately towards 
each other. We account for this in the model by stretching the 
right half of the King profile (long dashed line) by a factor /, 
which has to be determined, i.e. the original profile /CKing(S) is 
replaced with A-stretch(0) = '<'King(/^), < / < 1, for positive val- 
ues of 6. We will call / the "stretch factor". The contribution of 
the filament (dotted line) /Cfii(0) is added to the stretched King 
profile. The result is the observed surface mass density profile 
on the right-hand side (short dashed) which can be described 
by 
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Fig. 18. 12*^*1 maps of toy models of two clusters. Top panel: without 
a connecting filament. Bottom panel: with a filament running between 
both clusters. Quadrupole moments in the left panel were computed 
in a 2' aperture, in the right panel in a 3' aperture. The radii of the 
apertures correspond to the radii of the circles at the lower left corners 
of the respective panel. White contours are surface mass density; the 
black circles correspond the ^200 of the clusters. 

Since the mass profile of the filament alone is not accessible 
by observations, we have to determine a point on the observed 
mass profile Ko^J^O) that we treat as the "end" of the cluster and 
the "start" of the filament. We tried this using the following 
procedure: The unstretched King profile, observed on the left- 
hand side where < 0, is stretched by the factor /, to model the 
influence of tidal stretching. By this step we try to obtain the 
(unobservable) cluster profile A'stretch(S) on the right side without 
the contribution of the filament. 

This stretched profile is then compared to the observed pro- 
file /Cobs containing the contribution of the filament by comput- 
ing the goodness of fit 



X 



=Z 



/ ^stretch(g<)-^obs(gi) 

\ a-: 



(7) 



Kohsiff) = '^King(/0) + Kf,\{0) . 



(6) 



at sample points 0, in the reconstruction along the main axis of 
the system. Typically, the spacing of the sample points will be 
that of the grid on which the reconstruction was performed. 
Linear interpolation between grid points will be used if the 
sampling points do not exactly coincide with the grid points. 
cr, is the estimated error in /Cobs at the /th point. We define our 
sample points such that Oq - 0, i.e. 60 is placed at the clus- 
ter center. Usually, we will set /q = 1 ■ Note that in our model 
the contribution from the point 9 - Q always vanishes as by 
definition the observed and the stretched profile have the same 
value. 

X^ is repeatedly computed for increasing values of A^. We 
can define the "end of the cluster" and the "start of the fila- 
ment" by the point 0cut = Sn, where the probability that /Cstretch 
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Fig. 19. Simple model of the surface mass density distribution of an 
elliptical cluster and a filamentary extension along the main axis of 
the system. The solid line is a symmetric King profile, the long dashed 
line is the same King profile stretched by a factor to introduce the 
ellipticity seen in simulation. The filament is modeled as a separate 
component (dotted line). The observed profile (short dashed line) is 
the sum of the filament and the stretched King profile. The axes are 
labeled in arbitrary units. The vertical lines exemplify typical values 
for 6(i, and O^^i; see text for details. 



is a good representation of /Cobs falls below a pre-defined level, 
which we call the "cut-off confidence level" xluv 

We now have to find a way to determine the stretch factor 
/. For this, we assume that a position 9^, < < Ofn < Ocut 
exists, such that the influence of Kf\i(0 < Ofn) is negligible, i.e. 
we assume that the observed profile is a fair representation of 
the (unobservable) stretched profile /fstretch- The stretch factor / 
can then be determined by fitting the unstretched profile, which 
we obtain from observations at 6 < Q, to the inner portion (i.e. 
at < Ofit) of the observed profile. This "stretch factor fit" was 
done using ax^ minimization. 

The "cut-oflF parameter" 0cut and the "cut-ofi" confidence 
level" ;if^m have to be determined from simulations. Fig. 20 
shows the mass profiles to the left and right of the center of 
the left cluster in the reconstruction displayed in Fig. 13 along 
the main axis of that system. For simplicity the error bars were 
assumed to be equal to the standard deviation of a reconstructed 
mass map of a randomized catalog of background galaxies. 

We determined several combinations of the cut-off param- 
eter Ocut and confidence level ;^^^m that match the visual impres- 
sion of filament beginning and cluster end. However, if these 
were applied to clusters from other simulations, the separation 
point between cluster and filament was placed at non-sensical 
positions. 

We also modified the starting position 0,„ in the summation 
in eq. (7). First, we placed it at Ofn in order to exclude the cen- 



Fig. 20. Surface mass density profiles of the cluster on the left in the 
reconstruction displayed in Fig. 13 along the line connecting both 
cluster centers. The crosses mark the surface mass density in the fila- 
ment part, the dashes the surface mass density on the lefthand side of 
the cluster. The x-axis denotes the distance from the cluster center in 
arcminutes. 



tral region, which by definition of this procedure has a small 
X^- Second, we calculated ;^^^ in a moving window of fixed size 
and set the separation point between cluster and filament to the 
start of the window for which x^ fell below the cut-off confi- 
dence level. This was done for various window sizes and confi- 
dence levels. Again, parameters that worked well for one clus- 
ter failed completely for others in both approaches. 

5.5. Quadrupole moment map of A 222/223 

Having found in the previous section that aperture quadrupole 
moments in principle can be used to quantify the presence of 
a filamentary structure, if we are able to choose the right size 
of the aperture, we now apply this method to the A 222/223 
system. Fig. 21 shows a \Q^^^\ map with a weight function with 
3f2 radius. The white significance contours show a quadrupole 
moment signal on the filament that reaches a peak SNR of 3.0. 
The filter scale was chosen to be the same in which the aperture 
mass statistics gave the most significant signal in the interclus- 
ter region. The I2*^'| signal does not fully trace the filament 
candidate but only the Western part of it and an extension to- 
wards the mass peak in the East. The peak SNR is most likely 
enhanced by the trough at the Western edge of the mass bridge. 
It is not surprising that the significance of the quadrupole mo- 
ments on the possible filament is not very high. The aperture 
mass statistics already gave a relatively low SNR. The AMM 
statistics uses data within the same aperture as M^p but gives 
more information, namely instead of the mass in the aperture. 
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Fig. 21.Aperture quadrupole moment map 
of A 222/223 in an aperture of 3'2 radius, 
corresponding to the radius of the circle at 
the lower left corner. The thick white lines 
are SNR contours for |2*^'|, the lowest con- 
tour being at 2 and higher contours increas- 
ing in steps of 0.5. The aperture quadrupole 
moment on the "filament" region reaches a 
peak SNR of 3.0. The black lines are the 
contours of the mass reconstruction. The 
dashed circles indicate the ^200 from our best 
fit NFW models. 



it gives the mass distribution. Being generated from the same 
information, this naturally comes with a lower SNR. 

Several other features - most of them associated with the 
slopes of the two massive clusters - are also seen in Fig. 21. 
Interesting in the context of quantifying filaments is the I2*^'| 
statistics in the region between A 222 and the newly detected 
cluster SE of it. Indeed we see a signal with a peak signal-to- 
noise ratio of 2.7 extending between the two clusters. The mass 
reconstruction in Fig. 2 also shows a connection between both 
clusters but at a level that is dominated by the noise of the k- 
map. The redshift difference between A 222 and the new clus- 
ter, inferred from the color-magnitude diagram (Fig. 7), makes 
it very unlikely that these two clusters are connected by a fil- 
ament. It is much more probable that we are in the situation 
depicted by Fig. 16 in which the influence of the two individ- 
ual clusters leads to a quadrupole moment in the intercluster 
region. 

6. Discussion and conclusions 

Based on observations made with WFI at the ESO/MPG 2.2 m 
telescope we find a clear lensing signal from the Abell clus- 
ters A 222 and A 223. Comparing our lensing analysis with the 
virial masses and X-ray luminosities, we find that A 222/223 
forms a very complex system. Mass estimates vary consider- 
ably depending on the method. Assuming the best-fit NFW 
profiles of Sect. 3.2 /W20o(A 222) = 3.0!°'^ x lO^^h-^^MQ 
and /W20o(A 223) = 5.3![-^ x lO'^^/z^o'/W©. The masses of the 
best-fit SIS models within r2oo as determined from the NFW 
fit are higher for both cluster but compatible within their re- 



spective error bars: /Wsis(A 222) = 4.8+}'^ x lO^^/z^^/W© and 
/Wsis(A 223) = 7.3+}-^ x lO^'^h-l^MQ. These mass estimates are 
considerably lower than those derived from the virial theorem 
for an SIS model. Using the velocity dispersions from D02 we 
find /Wvir(A222) = 9.7+[-^ x lO^^h-^MQ and /Wvir(A 223) = 
12.0!|^ X 10^*h-^MQ. 

The M/L ratios we found in Sect. 4 are lower than the 
ones determined by D02 of M/Lr = (202 + 43)/z7oMo/Lo 
for A 222 and M/Lr = (149 + 33) /zvo Mq/Lq within a radius 
of 1.4/z"' Mpc but agree within the error bars of our values 
for the M/L ratios determined from the SIS model masses, and 
in the case of A 223 also with the M/L ratio from the NFW 
model. Two competing effects are responsible for this differ- 
ence. First and foremost, the weak lensing masses are lower 
than the masses D02 used. Second, also the luminosities deter- 
mined are lower than in D02. This has two reasons. First, D02 
analyze the Schechter luminosity function; this allows them to 
estimate the fraction of the total luminosity they observe, while 
we limit our analysis to the actually observed luminosity. Also, 
D02 correct the area available to fainter objects by subtracting 
the area occupied by brighter galaxies, which might obscure 
fainter ones. Both differences mean that we probably underes- 
timate the total luminosity of the clusters. Our M/L ratios are 
already at the lower end of common M/L ratios. A higher lu- 
minosity would lead to even lower M/L value making A 222 
and A 223 unusually luminous clusters, considering their mass. 
Although this system is complex and probably still in the pro- 
cess of collapsing, the M/L ratios of both clusters are very sim- 
ilar and do not exhibit variations like those observed by G02 in 
A 901/902. Variations between mass, optical, and X-ray lumi- 
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nosity are seen on smaller scales in A 223. A 223-N is very 
weak in the X-ray image, while it is the dominant sub-clump 
in the mass and optical luminosity density map. The latter may, 
however, be affected by the color selection that misses many of 
the unusually blue bright galaxies in A 223, especially in the 
Southern sub-clump. 

The weak lensing mass determination depends on the red- 
shift of the FBG which we assumed to be zfbg = 1 ■ This as- 
sumption is based on the redshift distribution of the Fontana 
et al. (1999) HDF-S photometric redshift catalog. Changes 
in the redshift distribution could change the absolute mass 
scale while leaving the dimensionless surface mass density 
and hence also the significance of the weak lensing signal un- 
changed. However, the A 222/223 clusters are at comparably 
low redshift and changes in zfbg affect the mass scale only 
weakly. To bring Msis(A 222) to the value of Mvir(A 222), the 
mean redshift of the faint background galaxies would have to 
move to Zfbg = 0.3. This is clearly unrealistic given the depth 
of our WFI images and the color selection we made. Although 
we cannot exclude deviations from the redshift distribution of 
Fontana et al. (1999), it is much more realistic to attribute the 
dififerences between virial and weak lensing masses to intrinsic 
cluster properties. 

From the visual impression of the galaxy distribution it is 
already obvious that this system is far from being relaxed. This 
can aff^ect the measured masses in several ways: First, the de- 
viation from circular symmetry certainly implies that the line 
of sight velocity dispersion is not equal to the velocity disper- 
sion along other axes in the clusters. If the clusters are oblate 
ellipsoid with their major axis lying along the line of sight, the 
measured velocity dispersions will overestimate the velocity 
dispersion. Second, if the clusters are not virialized, estimating 
their masses from the virial theorem of course can give signif- 
icant deviations from their actual mass. Finally, we could only 
successfully obtain weak lensing mass estimates with spheri- 
cal models, which are probably not a good representation of 
the actual system. Although both clusters are clearly ellipti- 
cal, fits with SIE models could not reliably reproduce the ob- 
served cluster properties. This does not come as a total sur- 
prise; King et al. (2002) already noticed that the shear log- 
likelihood function is much more sensitive to changes in the 
slope than to a possible cluster ellipticity. The insensitivity of 
the log-likelihood function to the ellipticity parameters means 
that the fitting procedure rather changes other cluster parame- 
ters than reproducing the actual ellipticity which we see in the 
parameter-free weak lensing reconstruction. This behavior and 
the difficulty to accurately fit elliptical models to shear data is 
confirmed by our simulations in Sect. 5.2. Although our as- 
sumption about the dispersion of intrinsic galaxy ellipticities 
and number density were much more optimistic than justified 
by our data, we could not recover the ellipticity and orientation 
of the clusters in the A'-body simulation. 

We found that the concentration parameter c of the NFW 
profile is poorly constrained if we omit the central regions of 
the cluster in order to avoid contamination with cluster galax- 
ies. This does not significantly affect the masses determined 
from fitting NFW profiles and was not of prime importance to 
the work presented here. From varying the radius of the circles 



in which shear information was ignored, we saw that galax- 
ies closer to the cluster center constrain the concentration pa- 
rameter better than those at large distances from the clusters. If 
one wants to determine concentration parameters more reliably, 
the fitting procedure has to be extended to include background 
galaxies close in projection to the cluster centers, while ensur- 
ing that faint cluster galaxies do not have a strong influence on 
the shear signal. 

The lensing reconstruction shows a "bridge" extending be- 
tween both clusters of the double cluster system. We devoted 
much effort to developing a method that could objectively de- 
cide whether this tantalizing evidence is indeed caused by a 
filament like it is predicted from A^-body simulations of struc- 
ture formation. Unfortunately, this was mostly done without 
success. The aperture quadrupole moment statistics in princi- 
ple has the power to detect the presence of a filament-shaped 
structure. To objectively apply it, one however needs to be able 
to separate clusters from the filaments connecting them. We did 
not find an objective way to do this and had to resort to subjec- 
tively defining the sizes of the apertures used. 

We would like to stress that this is not a problem of the 
weak lensing technique but stems from the fact that the descrip- 
tion of the cosmic web as filaments and galaxy clusters is based 
on the visual impression of A^-body simulations. Attempts to 
objectively separate these two components from each other re- 
quire a mathematical description which we tried to develop in 
Sect. 5.4. This was mostly unsuccessful because we could not 
find a procedure that reliably reproduces our visual impression. 
The visual impression of what a filament is, is often sufficient 
in simulations or redshift surveys where filaments stretching 
long distances between clusters are seen. In the case of close 
pairs of clusters - where we can hope to see filaments with to- 
day's telescopes - a more objective criterion is important, but 
difficult to find. 

We have not addressed the question how to distinguish the 
aperture quadrupole moment of a filamentary structure from 
that of a pure double cluster system in Sect. 5.3. We found 
that the quadrupole moment in a system with a filament ex- 
ceeds that of a halo-halo system without filament. Closer in- 
spection reveals that the shape of the quadrupole moment in the 
intercluster region changes if a filament is added to a two halo 
system. Because one can compute significances for AMM in a 
limited spatial region, a significant deviation from the expected 
shape of the quadrupole moment from a pure halo system could 
possibly be used to overcome this difficulty. This can only work 
if the signal-to-noise ratio of the aperture quadrupole moment 
is high. Possibly stacking several cluster pairs could provide a 
sufficiently high SNR. This could in principle be tested with 
our A^-body simulations but is beyond the scope of this paper 
in which we try to develop a criterion to quantify the evidence 
for filaments in single systems, like the A 222/223 system at 
hand. 

What can we then say about a possible filament between 
A 222/223? All observations presented in this work - weak 
lensing, optical, and X-ray - show evidence for a "filament" 
between the two clusters. The most compelling evidence prob- 
ably comes from the number density of color-selected early- 
type galaxies, which is present at the Icr level (Fig. 9). The 
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spectroscopic work of D02 and Proust et al. (2000) confirmed 
the presence of at least some galaxies at the cluster redshift in 
the intercluster region. Obtaining a larger spectroscopic sample 
in the intercluster region would allow us to spectroscopically 
confirm the significance of this overdensity and could provide 
insights into the correlation of star formation rates and matter 
density (e.g. Gray et al. 2004). The X-ray emission between the 
clusters is aligned with the overdensity in galaxy number and 
luminosity density. This provides further evidence for a "fila- 
ment" extending between A 222 and A 223. 

The signal level of the possible filament in the weak lensing 
map Fig. 2 is rather low compared to the clusters. The aperture 
quadrupole statistics has a signal at the 3cr level on the filament 
candidate but this signal may already be contaminated by the 
outskirts of the cluster in the aperture. The most striking "fea- 
ture" of the mass bridge seen in the k map is the misalignment 
with respect to the possible filament seen in the optical and X- 
ray maps. This can be interpreted in several ways. It could sug- 
gest that the surface mass density on the "true filament" defined 
by the position of the optical overdensity and X-ray emission 
is below our detection limit and what we see in the k map is 
a noise artifact. This possibility aside, the observed misalign- 
ment can have several causes. First, as we already discussed in 
Sect. 4, the influence of the many reflection features around the 
bright star West of A 223 on the weak lensing reconstruction is 
difficult to determine. It seems that the cluster peaks are shifted 
preferentially away from the reflection rings. The same could 
be true for the "filament" in the reconstruction. Second, the po- 
sition of structures inferred from weak lensing is affected by 
the noise of the reconstruction. This is especially true for low 
mass structures and is illustrated by our simulations using SIS 
models to infer the positional uncertainty of weak lensing re- 
constructed peaks in Sect. 3.2. It is possible that at least part of 
the observed misalignment is caused by the noise of the weak 
lensing method. Finally, one could in principle imagine that the 
off-set is real and a misalignment of dark and luminous mat- 
ter is present. This would require complex and possibly exotic 
physical processes that cause galaxies to form next to a dark 
matter filament and not in it. At present there is no good obser- 
vational support for such a scenario. We should note, however, 
that a misalignment between mass and light is also present in 
the filament candidate of G02. 

As we have not found an objective way to define what a 
filament in a close double cluster pair is, the question whether 
what we observe in A 222/223 constitutes a filament or not can 
also not be answered objectively. Thus, our filament candidate 
is - in this respect - not very different from those of Kaiser 
et al. (1998) and G02. The projected virial radii of the clus- 
ters marginally overlap. However, (1) the redshift difference 
between the clusters make an actual overlap of the clusters un- 
likely; (2) the projected mass in clusters falls off steeply, and 
weak lensing is currently not capable of mapping the cluster 
mass distribution out to the virial radius. A signature of a fila- 
ment should thus already be present inside the virial radius. 

The unambiguous weak lensing detection of a filament be- 
tween two clusters would provide a powerful support for the 
theory of structure formation and the "cosmic web". Taking the 
3cr signal of the quadrupole statistics on the filament candidate 
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Fig. A.l. Arclet candidates around the cD galaxy in A 222. North is 
up and East is to the left. The scale is 10" long. 



at face value, an increase of the number of density of FBGs by a 
factor of 2.8 could give a 5cr detection. Such number densities 
can be reached by 8 m class telescopes. The A 222/223 system, 
being only the third known candidate system to host a filament 
connecting two cluster, would be a good target for such a study. 
In fact, a weak lensing study of A 222/223 using SuprimeCam 
at the Subaru telescope is already underway (Miyazaki et al., 
in preparation). 

In addition to the lensing signal from Abell 222 and 
Abell 223 we found a significant mass peak SB of A 222. 
This peak coincides with an overdensity of galaxies. The color- 
magnitude diagram of these galaxies suggest that this newly 
found cluster is at a redshift of z ~ 0.4, but this estimate comes 
with a considerable uncertainty and requires spectroscopic con- 
firmation. A maximum likelihood fit to the shear data around 
this mass peak leads to a best-fit SIS model with a veloc- 
ity dispersion of 728;^[2o kms"'. This serendipitous detection 
again illustrates the power of weak lensing as a tool for cluster 
searches. 
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Appendix A: Strong lensing features in A 222 

Already in 1991 Small et al. (1991, SEF) found two candi- 
date arclets in the center of A 222. We also see two possible 
arclets in the center of A 222 displayed in Fig. A.l. Arclet 1 is 
the same as found by SEF and labeled A 222-1. Unfortunately, 
SEF's second candidate is not marked on the plate in their pa- 
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per, and as SEF give only distances from the cluster center and 
no position angle we do not know whether their second can- 
didate corresponds to ours. A comparison of the arclet candi- 
date properties between SEF and our candidates is given in ta- 
ble A.l. The distance measurements for A 222-1 and arclet 1 

Table A.l. Arclet candidate properties from SEF and our data. The 
column entries are distance from the center of the cD galaxy, axis 
ratio, and position angle measured clockwise from the north direction 
in Fig. A. 1 . 



Arc ID 


d/arcsec 


b/a 


pos. angle 


SEF 










A 222-1 




12.2 


4.7 


- 


A 222-2 




14.1 


2.6 


- 


this work 










arclet 1 




12.7 


2.8 


42° 


arclet 2 




10.8 


3.0 


140° 



are in good agreement but the values for the axis ratio show a 
clear deviation. The difference may be due to the comparably 
poor image quality in the work of SEF and blending with the 
nearby object to the South- West of arclet 1 . However, it must 
also be mentioned that the determination of the axial ratio is 
relatively uncertain and we estimate its error to be of the order 
^0.6. 

Given the discrepancy between the distance measurements 
for A 222-2 and arclet 2 it is unlikely that these are the same 
objects. 

Unfortunately, the V band image is not deep enough to 
show the candidate arclets, so that no color information is avail- 
able. 

Appendix B: Multipole moments 

Schneider & Bartelmann (1997) define the complex nth-order 
aperture multipole moment as 



gW(6lo)= d'-eff'U(\0\)e""f 400 + 0) 



(B.l) 



with a radially symmetric weight function L'^d^l). For n = 2 
eq. (B.l) expresses the aperture quadrupole moment in terms 
of the surface mass density. Based on this definition, an ex- 
pression for the aperture moments in terms of shear estimates 
may be found (Schneider & Bartelmann 1997): 



1 ^ 

n ^-^ 

1=1 



eflnUiOi) + 0U'(ed] 



(B.2) 



where n is the number density of galaxies in the aperture, 
(6i, (pi) are the polar coordinates of the /th galaxy with respect to 
^0, and Sti - -'Rieit'^^^') and Sxi - --^(e,e"^"^0 are the tangen- 
tial and cross components of the shear estimate, respectively, 
with respect to ^o- Here U'(ff) is the derivative of the weight 
function. 



We now show that the definition (B.l) cannot be general- 
ized to non-radially symmetric filters U{0). We partially inte- 
grate eq. (B.l) with respect to if and obtain 
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where for simplicity we have set 00-0 without loss of gen- 
erality. The integral over the first term in this expression can 
be expressed in terms of the shear in analogy to Schneider & 
Bartelmann (1997), while we integrate the second term again 
by parts, this time with respect to 9. This integration over the 
second term leads to an expression similar to (B.3) with inte- 
grals over two terms; one that can be readily expressed in terms 
of the shear, the other requiring further integration by parts, 
and so on. Eventually, the aperture multipole moment can be 
expressed in terms of the shear as an infinite series of integrals: 
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jr^\n/ Jo Jo oOdip' 



with the requirement on the weight function that 



rx0"' 



d_ 
dip 



U(0) -» 0, for -^ and ^ oo 



(B.5) 



so that the integrals exist, yt and jx are the tangential and cross 
components of the shear, in analogy to St and Sx above. It turns 
out that this sum in general does not converge. E.g. for n — 2 
and 



U{0) = f{0) 



'*!(• 



2i^ 



-He 



-2i^l 



(B.6) 



where f{6) is a smooth, positive, and finite function satisfying 
the condition (B.5), the sum in eq. (B.4) oscillates around 0. 
This behavior can be understood if we insert (B.6) into (B.l) 
and integrate by parts, while setting k{0) to a constant value ko- 
The quadrupole moment then depends on /cq, unless we allow 
9^f{9) to be compensated. A constant ko does not influence the 
shear. Hence, the sum in eq. (B.4) cannot converge. We thus 
find that the mass-sheet degeneracy prevents us from comput- 
ing aperture multipole moments in non-circular apertures. 
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